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1 .  INTRODUCTION 

The  following  interim  technical  report  has  been  prepared  by 
SAD  Dynamics,  Inc.  under  Contract  No.  DAAG29-83-C-0004  to  the 
U.S.  Army  Research  Office,  Research  Triangle  Park,  NC. 

This  interim  report  presents  the  detailed  formulation  of  the 
equations  of  motion  of  a  projectile  of  finite  geometric  and  iner¬ 
tial  properties  traveling  in  a  flexible  gun  tube.  The  formulation 
permits  the  projectile  six  degrees-of-freedom  relative  to  the  gun 
tube;  three  orthogonal  translational  motions  of  the  projectile 
c.g.  relative  to  the  instantaneous  gun  tube  axis  and  three 
(Eulerian)  rotations  of  the  projectile  about  its  c.g.  (related 
to  projectile  pitch,  yaw  and  roll  motions).  The  formulation 
accounts  for  projectile  spin,  mass  eccentricity,  elastic  deforma¬ 
tion  of  the  projectile  rotating  band  and  bourrelet,  and  projec¬ 
tile/bore  interfacial  friction  and  torque  transmission.  Further¬ 
more,  the  formulation  is  unrestricted  regarding  the  nature  of  gun 
tube  motion. 

The  detailed  formulation  is  presented  in  Section  2.  A 
comparison  of  this  formulation  with  other  projectile  descriptions 
and  formulations  is  presented  in  Section  3»  Conclusions  are 
presented  in  Section  k.  References  cited  are  presented  in 
Section  5.  The  Euler  angles  and  associated  transformations 
employed  are  defined  in  the  Appendix. 
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2.  FORMULATION  OF  PROJECTILE  EQUATIONS  OF  MOTION 

The  simultaneous  system  of  differential  equations  which 

describe  the  general,  six  degree-of-freedom  motion  of  a  projoctile 

\ 

traveling  in  a  flexible  gun  tube  are  formulated  in  this  section. 

The  pertinent  equations  are  formulated  in  terms  of  the  loads  and 

moments  applied  to  the  projectile  in  Section  2.1.  Simplification 

of  the  projectile  angular  velocity  and  acceleration  expressions 

is  discussed  in  Section  2.2.  The  applied  loads  and  moments  acting 

on  the  projectile  are  defined  in  Section  2.3*  The  solution 

technique,  which  allows  simultaneous  solution  of  the  equations 

herein  developed  with  the  equations  of  the  gun  dynamics  simulation 

code  previously  developed,  is  discussed  briefly  in  Section  2.4 

and  will  be  presented  in  greater  detail  in  a  separate  report. 

2. 1  Projectile  Equations  of  Motion  in  Terms  of  Applied 
Loads  and  Moments 

The  projectile  is  described  subject  to  the  following 
assumptions : 

(i)  the  projectile  is  assumed  to  consist  of  a  main  body, 
rotating  band  and  bourrelet; 

(ii)  the  projectile  main-body  is  assumed  to  behave  as  a 
rigid  body  of  finite  geometry  and  inertia; 

(iii)  the  rotating  band  and  bourrelet  are  assumed  to  behave 
elastically  (although  not  necessarily  linearly 
elastic ) ; 

(iv)  the  projectile  c.g.  is  assumed  to  be  eccentrically 
located  relative  to  its  geometric  center; 

(v)  the  projectile  is  permitted  six  degrees-of-freedom 
relative  to  the  gun  tube,  namely  pitch,  yaw  and  roll 
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about  Eulerian  axes,  and  translational  motions  along 
body-fixed  (orthogonal)  axes. 

The  following  coordinate  systems  are  introduced  for  the 
purpose  of  formulating  the  equations  of  motion  of  the  projectile: 

(i)  S  ,  with  coordinates  (x  ,  y  ,  z  )  and  unit  triad 

(£  ,  3*  ),  is  defined  as  an  inertial  reference 

frame  fixed  at  the  initial  projectile  position  such 
that  its  orientation  coincides  with  the  global 
(earth-fixed)  coordinate  system  defined  in  tie  gun 
dynamics  simulation  code; 

(ii)  Sq,  with  coordinates  (x0»yo*20)  and  unit  triad 

(£  ,  j*  ),  is  defined  as  an  intermediate  reference 
o  o  o 

frame  whose  origin  translates  with  the  projectile 
along  the  gun  tube  axis  and  rotates  with  the  projectile 
about  the  instantaneous  tangent  to  the  gun  tube  axis. 
iQ  is  directed  toward  the  gun  tube  muzzle,  along  the 
instantaneous  gun  tube  axis.  At  the  initial  projectile 
position,  3“0  is  directed  along  the  radius  from  the  gun 
tube  centerline  to  the  projectile  c.g.j 

(iii)  S,  with  coordinates  (x,y,z)  and  unit  triad  (t.j.fi), 
is  defined  as  a  body-fixed  coordinate  system  whose 
origin  is  located  at  the  projectile  c.g.  £  is 
directed  toward  the  projectile  nose,  parallel  to  its 
geometric  axis.  At  the  initial  projectile  position, 

J  is  parallel  to  J  . 

i 

It  is  noted  that  the  kinematic  relations  between  S  and  S 

o 

incorporate  gun  tube  motion,  as  well  as  the  two  degrees-of-freedom 
of  the  projectile  corresponding  to  translational  motion  along  the 


gun  tube  axis  and  spin  (both  of  which  are  prescribed  by  interior 
ballistics  data).  The  remaining  four  degrees-of-f reedorn  of  the 


projectile  are  formulated  in  S  relative  to  S  .  Referring  to 

Figure  1 ,  these  are  defined  as  the  translational  displacements 

of  the  projectile  c.g.  in  the  J  -I*  plane  (denoted  as  y  and 

oo  '  eg 

z  ),  and  the  Euler  angles  corresponding  to  projectile  pitch  and 
eg 

yaw,  namely  y>  and  9  respectively. 

For  the  purpose  of  formulating  the  velocity  and  acceleration 
of  the  origin  of  Sq  relative  to  S  ,  the  Euler  angles  (' f'0*90*(P0 ) 
defined  in  the  Appendix  are  introduced  subject  to  the  transforma¬ 
tion 

V-  K  l  -  ll  f.  *  £s  l  ] 

j'*  An  L  +  Allfi*  Au  ie  > 

%-  Alts  +  A\X  , 

where  (referring  to  the  Appendix), 

/„  =  cos  t.  cos  9. 

/*  »  cost,  sine,  Sinf.  -  Siiit,  cos  f, 

A\  -  cost,  sin.9.tosf.  +  sinyslnf. 

H,=  sinf.  cos  0, 

Ah*  sinyssine.sinf.+cosyco’sf, 
jtU-sintsine.cos'f.-cosr.sinf, 

A‘„*  -  sin  9, 

Ah  =  cosd.sin  9. 

/•„=cos6'Cosf. 
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Since  S  is  constrained  to  move  wLth  the  projectile,  the 
o 

position  vector  of  its  origin,  Rq,  as  shown  Ln  figure  2,  has 
components  which  are  functions  of  time,  t,  as  well  as  distance 
along  the  gun  tube  axis,  s.  Hence,  the  velocity  of  the  origin 
of  So  relative  to  S  ,  namely  vq  (=  dRQ/dt),  is  given  as 


if*  it  I ift  ™ 


Since  i  lies  along  the  instantaneous  tangent  to  the  gun 
o 

tube  axis  it  follows  that 

0°  M  f°  0* 

Xu 


00 


Hence , 


ft  -  [%*£■%]  *'♦  L  jO ] %  (5) 


Applying  the  transformation  presented  in  equation  (l),  and 
noting  that  ds/dt  denotes  the  instantaneous  projectile  velocity, 
Vp  (as  prescribed  by  interior  ballistics  data),  there  results 
from  equation  (5)  the  velocity  of  Sq  in  the  form 


Since  S  is  earth-fixed,  the  instantaneous  acceleration  of 
Sq  is  readily  obtained  by  differentiating  equation  (3)  with 
respect  to  time  noting  that 


ji. 


^  as 


(7) 


from  which  there  results 


6 


where  a  denotes  dv  /dt. 

P  P 

Once  again  applying  the  t  ransf ormation  presented  in  equation 

(l),  and  noting  equation  (M,  there  results  the  acceleration  of 

S  in  the  form 
o 


Noting  that  Sq  is  constrained  to  move  with  the  gun  tube  in 
the  Jo4i  plane,  it  has  angular  velocity,  U>0»  relative  to  S  .  The 
expression  for  (written  relative  to  Sq),  is  given  as 

CJ, .  *  l»  +  4  (1  1  ^ 
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where  (referring  to  the  Appendix), 


u*.  =  f.-VisinB.  ) 

ut'  =  6.  cos  f.  +  ycosGl,  sin  f.  * 

Qlt  =  -d.sinf.+  y',tosB.co%f. 

•  •  • 

and  %<B0  and  are  obtained  applying  equation  (7)* 
Noting  further  that 


Mj 

cLt 


(*)  +  to  X  (  ) 


t 

there  results  the  angular  acceleration  of  relative  to  S 


(13) 


(14) 


where  (from  equation  ( 12)), 


y-Vsine'-itbcose. 

Cj^  0oco$f.+f,co$6.sin<f.-i'emsin6.sinf.  +  ul0<%  > 

GJ^-GLsihf.+ycos&oCOsH  -  \6msih6^cosf0  - 


(15) 


Equations  (6),  (9)i  ( 1 1 )  and  (l4)  prescribe  the  linear  and 

angular  velocities  and  accelerations  of  SQ  relative  to  S  .  It 

remains  to  prescribe  the  motion  of  the  projectile  c.g.  relative 

to  S  . 
o 

Letting  v  and  a  respectively  denote  the  translational 
eg  eg 

velocity  and  acceleration  of  S  relative  to  S0#  it  follows  from 
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Figure  1  that 


V1  *  %f’  *  *'3^' 

S=  %f.  *  z.3-&  , 


(16) 


Letting  denote  tho  nn^lar  velocity  of  S  relative  Sq, 


and  further  specifying  that 


(17) 


W.J  =  W«,,  z  +  *  6J*.,-* 

the  components  of  cO  are  prescribed  (analogous  to  equation 

eg 

(12))  as  functions  of  the  Euler  angles  (Y',0,f>)  of  S  relative  to 


S  in  the  form 
o 


gj,cj=  f  - 

6c^tj  *  0  c<?3  f  ♦  y'cosOsin^P 
ojZc^  -0sin<P+V'co3dcos>y 


(18) 


Noting  that  Sq  is  prescribed  to  rotate  with  S,  equation 


(17)  must  be  subjected  to  the  constraint 


<£y  t.*o 


(19) 


Introducing  the  transf ormation  from  S  to  SQ,  namely 

l  -  Ai  *  4  As  £ 

=  AxX  *  +  A*  *L  ^ 

~  Ail *  +  An  £ 


(20) 


and  noting  that  the  direction  cosines,  , ,  are  defined  in  terms 
of  the  Euler  angles  analogous  to  the  definition  presented 

in  equation  (2),  there  results  from  equations  ( 1 7 )  thru  (20) 

Q  (21) 

which  prescribes  the  dependence  between  the  Euler  angles  of  S 
required  to  satisfy  the  constraint  presented  in  equation  (19). 

Equations  (l6)  and  (l7)  (subject  to  equation  (2l))  respec¬ 
tively  prescribe  the  translational  and  angular  motions  of  the 
projectile  c.g.  relative  to  Sq. 

Referring  to  Figures  1  and  2,  tho  instantaneous  position 

t  _ 

vector  of  S  relative  to  S  is  given  by  R  +r  .  Hence,  noting 

o  eg 

equation  (13),  the  (absolute)  translational  acceleration  of  the 
projectile  c.g.  (relative  to  S  ),  namely  a,  is  given  by 

a  =  at+ + ti'XKi  +  oj'XtQ'Xf:,)  (22) 

while  the  (total)  angular  velocity  of  the  projectile  (relative  to 

i 

S  )  is  given  by 

U)  =  U)0  +  (23) 

Applying  equation  ( 13)  to  each  vector  on  the  right-hand  side 
of  equation  (23)  there  results  the  (total)  angular  accelex'ation 
of  the  projectile  (relative  to  S  ) 

i)  =  CD0  +  CJt j  +  GJ0X  (24) 

Substituting  equations  (9),  (ll),  (l4),  and  ( 1 6 )  into  equa¬ 
tion  (22)  there  results 
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a.  =■  [a,.  +  V*  w'.y<»+  ^(^H.^-^.,’  <^*V  <*Wl  ’• + 

+  l-a>  *  V*) ' 2  “«•*•» "  V^--~  H.H,)*  < )]  £  + 

♦  [v  V*^V^+tK^*^‘ ■&  (25) 

Transforming  equation  (ll)  from  SQ  to  S  via  equation  (20) 
and  substituting  the  result,  along  with  equation  (l7)>  into 
equation  (23),  there  results 

CJ  =  GJ*  1  ■*■  CJj  ^  +  CJj 

where , 

CJ„  =  CJx,3  +  X  ♦  X,  +  X/ 

=  +  Xi  +  -Xi  +  X* 

CJ|  s  ^£tj  +  +  X»^;#  4  Xj^f, 

Transforming  equation  ( 1  ?« )  from  Sq  to  S  via  equation  (20), 

differentiating  the  components  of  equation  (17)  with  respect  to 

time,  and  substituting  these  results,  along  with  equation  (17), 

into  equation  (24),  there  results 

*  “J*  •  •  3"  <  . 

W  =  ^  t  +  +  tJ,#  (28) 

where , 
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U)%  *  in  ^  *4r,+  (  ^/J  +  ^»  Wjf,"+-^«l  W*J  +'| 

~  i-u. <J*o *  inty*  in. ) 

<V  <V  ^*<V  A<i  +A\ + J^/huJ  ♦  I  (29) 

~  ^C,J  (  in  ^*0  *  in  ^jf«  +A  ^#,  ) 

6J*  *  +  y^/j * in tyt+ in ty, *  (dftJititixS iiiWi'+in  "*" 

~~  ( i,x  +  i*i  *  in  )  j 

and  (from  equation  ( 1 8 ) ) , 

U)x^--  <f-y'sih9  -T9  cos9 

(j^  ecos9+Yco$esinf-y'6ym6siKi  +  Uiti<f  >  (30) 

kJetJ9  -  ds'in*?  +*tc osdcosf  ~y/9sindco$cP-  Uftjf 


and  (from  equation  (21 )) 


ip,  tsS^Q^-i^l 

•  cote  "  cost)  [ 


COS1)* 


flfanT'ta.n^]  (31) 


We  are  now  in  a  position  to  formulate  the  projectile  equa¬ 
tions  of  motion.  Letting  F  (which  remains  to  be  prescribed) 
denote  the  resultant  applied  force  acting  at  the  projectile  c.g. , 
and  applying  Newton's  law  of  motion  with  'a  as  prescribed  in  equa¬ 
tion  (25),  there  results  the  three  scalar  equations  of  transla- 


tional  motion  of  the  projectile  c.g.  (written  relative  to  SQ )  In 
the  form 


f^0=  +  w,," ^*#Wij#) -  (y,,-»£)(<J|,t+ 

£.*  »,[<V  2,+ZMt.y^  *  (y,,*£)tw..-*  u,  w,.)-  z«,  (<  ♦  w,*j] 


(32) 


where  m  denotes  the  mass  of  the  projectile,  and  P  ,  F  and  F 

P  X  V  25 

^  o  J  o  o 

denote  the  components  of  F  along  the  respective  axes  of  Sq. 

Letting  M  (which  also  remains  to  be  prescribed)  denote  the 

resultant  applied  moment  acting  about  the  projectile  c.g»,  and 

applying  the  principle  of  angular  momentum  with  Q  and  u)  prescribed 

respectively  in  equations  (2<>)  and  (28),  there  results  the  three 

scalar  equations  which  prescribe  angular  motion  of  the  projectile 

(written  relative  to  S  in  order  to  preclude  introducing  time 

derivatives  of  the  projectile  inertia  tensor)  in  the  form 

Mu  “  X^Uiy  OJt  ) 

M,  =  I,,^ *(I„-I„)w,H,+ I,.(w,<<y  w« r  (33) 

M,=  I„Wj )  , 

where  M^,  M  and  Mz  denote  the  components  of  M  about  the  respec¬ 
tive  axes  of  S,  and  I  ,  I  _,  ...»  I  denote  the  elements  of  the 

xx  yj  zy 

projectile  inertia  tensor  written  relative  to  S. 

Within  the  framework  of  the  assumptions  introduced,  equations 


(32)  and  (33)  present  the  gonorul  equations  of  j>roj«*i  l.lte  motion 


in  terms  of  the  applied  loads  and  moments • 

2 . 2  Simplification  of  Angular  Velocity  and 
Acceleration  Expressions 

The  expressions  for  the  terms  representing  the  angular 
velocity  and  acceleration  components  entering  equation  (33)  may 
be  greatly  simplified  by  noting  that  for  most  practical  applica¬ 
tions  f and  6  aro  sufficiently  small  such  that 


sin'f  *  y  * 


cos V 


SinO  a*  Q  ;  cos  6^1  J 


(3*0 


Under  this  condition  it  follows  from  equation  (21 )  that 

<f  ^  'y  e  (35) 

while  from  the  mean  value  theorem  it  follows  that  the  integral  of 
of  equation  (35)  will  at  most  be  of  the  order  Y9 .  Hence,  to 
first  order 

y  0  (36) 


Imposing  the  above  conditions  and  retaining  only  linear  (first 
order)  terms  inland  0,  the  direction  cosines  for  the  transforma¬ 
tion  from  S  to  SQ ,  as  given  in  equation  (20),  simplify  to 


(37) 
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Noting  that  equations  (l8)  and  (30)  simplify  respectively  to 


**  td-Gy  1 


Ct^€j  ~  Q 

&  i' 


and 


60y  ^  >  0  “  6  V* 

I 

tk  ^  e  *  yyo 

**3 

ut,^  y' -  eey- -y-e1 


(38) 


(39) 


and  substituting  these  expressions,  along  with  equation  (37)  into 
equations  (27)  and  (29),  there  results  the  greatly  simplified 
expressions  for  the  angular  veloc i ty  and  acceleration  components 
entering  equation  (33)  in  the  form 


cj%  »  u)*t  ♦  y  ( u/]  +  6 )  -  d(Ut  +  y ) 

=  U)^  t  9  *  y  OJ*. 


(40) 


and 


16 


4-  +  ' 

ty* £))  +  e +v($V'-u%m)-y'(<4)*sY'i*))')+y'0(ecju+u>u)  /  (**i) 


2.3  Applied  Loads  and  Momenta 

The  equations  of  projectile  motion  previously  developed 
require  specification  of  the  applied  loads,  F  (written  relative 
to  ■SQ),  and  the  applied  moments,  M  (written  relative  to  S),  for 
completion  of  the  formulation.  These  loads  and  moments  arise  as 
a  consequence  of  the  projectile  weight,  interfacial  contact  of 
the  rotating  band  and  bourrelet  with  the  bore,  propellant  gas 
pressure  acting  at  the  base  of  the  projectile,  and  compressed 
air  ahead  of  the  projectile. 

2.3.1  Projectile  Weight  Loading 
The  load  applied  to  the  projectile  c.g.  due  to  its  own 
weight  is  given  by 


where  m^  and  g  denote  respectively  the  projectile  mass  and 
gravitational  acceleration. 

Applying  the  transformation  given  in  equation  (l),  equation 
(42)  is  written  in  component  form  (relative  to  SQ)  for  applica¬ 
tion  to  equation  (32),  as  follows 
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rw= 

r> 

r"W 

t  3 


s  -  5/n  co5  0* 

-Mpjlsiny  sinG'S'intf*  cos^  cos?,)  > 
->hp<)($in'%  SinGLcostf  -  cos^sinU) 


(43) 


2.3*2  Rotating-Band/Bore  and  Bourrelet/Boi'e 
Interfacial  Contact  Loadings 

The  loads  and  moments  applied  at  the  projectile  c.g.  due  to 
interfacial  contact  of  the  rotating  band  and  bourrelet  with  the 
bore  are  derived  subject  to  the  following  assumptions} 

(i)  radial  deformation  of  the  rotating  band  and  bourrelet 
are  characterized  by  Winkler  foundation  models; 
rendering  rotating-band/bore  and  bourrelet/bore 
interfacial  load  distributions  which  are  directed 
radially,  with  local  magnitude  determined  by 
resultant  local  radial  displacement; 

(il)  both  the  rotating  band  and  bourrelet  are  permitted 
radial  pre-loaded  initial  states; 

(iii)  the  instantaneous  annular  segments  defining  rotating- 
band/bore  and  bourrelet/bore  intorfacial  contact  are 
kinematically  determined  by  the  extent  of  annular 
compression  within  the  rotating  band  and  bourrelet 

at  each  instant; 

(iv)  the  projectile  main-body  is  assumed  to  be  rigid 
compared  to  the  rotating  band  arid  bourrelet; 

(v)  curvature  of  the  gun  tube  axis  between  the  planes 
containing  the  rotating  band  and  bourrelet  is 
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neglected; 

(vl)  rifling  torque  is  transmitted  without  slippage  via 
a  uniformly  distributed  circumferential  load  acting 
over  the  instantaneous  annular  segment  of  the 
rotating  band  in  contact  with  the  bore; 

(vii)  the  Euler  angles  y  and  Q  satisfy  equation  (34). 

In  view  of  Assumptions  (iv)  and  (vii),  the  displacement 
relative  to  Sq,  cf,  of  any  point  of  the  projectile  in  a  plane 
perpendicular  to  its  geometric  axis  is  given  by 

2  =  u.,-JL9)i.  (w) 

where  X  denotes  the  perpendicular  distance  from  the  projectile 
c.g.  to  the  plane  of  interest.  It  is  noted  that  Jt>0  implies 
that  the  plane  of  interest  is  forward  of  the  projectile  c.g., 

JL  m  0  implies  'that  the  plane  of  interest  contains  the  projectile 
c.g.,  while  l  0  implies  that  the  plane  of  interest  is  aft  of  the 
c.g. 

The  radial  component  of  this  displacement,  is  given  as 

where  denotes  the  angle  between  the  projection  of  the  y^-axis 
onto  the  plane  of  interest  and  the  line  from  the  gun  tube 
centerline  to  the  point  under  consideration. 

The  maximum  radial  displacement,  <£max»  in  the  plane  of 
interest  is  determined  by  setting 

=  0 
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Letting  *P  denote  the  orientation  of  maxim  am  radial 
a  *  max 

displacement ,  there  results  from  equation  (45) 


and  hence, 


<Z>  „  J.  -/  /  “  /■  &  ) 

ran  (  ^  +  / 

<£.„=  /fy,, 


(**7) 


(48) 


In  view  of  equations  (47)  and  (^8),  there  results  from 
equation  (45) 


Sr  =  cL,  cos(?-1L,) 


(«  9) 


Equation  (49)  denotes  the  radial  component  of  displacement 
of  any  point  in  the  plane  of  interest,  and  in  particular,  it  may 
be  used  to  specify  the  radial  displacement  of  the  planes  containing 
the  rotating  band  and  bourrelet  at  the  boro  surface.  Tt  should 
be  noted  that  the  radial  displacement  as  given  by  equation  (49) 
is  symmetric  with  respect  to  the  angle  lPmax 

In  the  sequel,  quantities  referred  to  the  plane  containing 
the  bourrelet  are  assigned  the  subscript  " 1 M ,  quantities  referred 
to  the  plane  containing  the  rotating  band  are  assigned  the 
subscript  "2".  Furthermore,  ^lnax  is  replaced  by  either  or 
and  tfmax  is  replaced  by  either  f,  or  f2. 

Assuming  that  the  bourrelet  is  forward  of  the  projectile 
c.g.,  while  the  rotating  band  is  aft  of  the  c.g.,  the  radial 
displacement  at  the  bourrelet  due  to  projectile  motion  is  given 
by 

Jr(  =  Slcos(f-'fl  )  (50) 
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where , 


(51) 


while  the  radial  displacement  at  the  rotating  band  is  given  by 


L  =  -£) 

Mf 


(52) 


where  , 


4  ( z,t 

<J>  -  4-  -I  /  ^  ) 

7i-  tan  J 


(53) 


Letting  and  <5°  respectively  denote  the  initial  radial 
pre-compressions  of  the  bourrelet  and  rotating  band  when  the  bore 
centerline  and  projectile  geometric  axis  are  initially  aligned, 
it  follows  from  Assumptions  (i)  thru  (iv)  that  the  radial  loads 
per  unit  circumferential  length  are  given  as 


Ri<f)  =  i6i[<J,'  + J  i  =  /,z 


(5M 


where  (1=1,2)  respectively  denote  the  radial  spring  stiffness 
per  unit  circumferential  length  of  the  bourrelet  and  rotating 
band. 

It  is  noted  that  the  expressions  presented  in  equation 
(54)  are  valid  for  compression  only,  that  is,  for 


+  COS('f-'fi)]-  0  ■,  i -1,2  (55) 
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Defining  f  (is  1,2)  as  the  instantaneous  maximum  annular 
Ci 

range  (measured  from  for  compressive  loading  of  the  bourrelet 

and  rotating  band,  equation  (5*0  may  be  written  in  the  form 

RiW--<  <56> 

l  0  i  (%-%)-? HVi+%) 

It  follows  from  equation  (55)  and  the  definition  of  <f  that 
o  Ci 

when  ^  =  0,  <pc  =  TT/ 2.  Hence,  in  the  absence  of  pre-compression, 

loading  takes  place  on  half  of  the  bore  surface,  and  this  loading 

is  symmetric  with  respect  to  the  angular  orientation  of  maximum 

radial  displacement.  When  £?  s  <£.  it  follows  that  =  TT.  Hence, 

li  '/ci 

when  displacement  due  to  projectile  motion  is  less  than  the  pre- 
compression,  contact  is  maintained  over  the  full  bore  surface. 

For  intermediate  cases,  namely  0  <  <  8 i ,  it  follows  that 

^  =  tan' 

The  radial  load  distribution  corresponding  to  an  intermediate 
case  is  depicted  in  Figure  3. 

The  instantaneous  interfacial,  contact  loads  acting  on  a 
differential  element  of  either  the  rotating  band  or  bourrelet  are 
as  depicted  in  Figure  4. 

Under  Assumption  (vi),  the  friction  force,  A  acts  tangen¬ 
tial  to  the  rifling  surface  in  the  direction  opposing  projectile 
motion  at  the  surface.  The  angle  ot  depicted  in  Figure  4  is  given 
by 

oC  =  ta.n  1  ( Ik)  (38) 
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Figure  3  - 


where  r,_  and  t*  denote  respectively  the  bore  radius  and  rifling 
b 

twist.  In  addition,  in  accordance  with  Assumption  (vi)  the 
torque,  T,  applied  to  the  rotating  band  per  unit  circumferential 
length  is  assumed  to  be  a  constant  independent  of  over  the 
range  of  compression. 

Letting  dF^  denote  the  incremental  resultant  contact  load 

acting  on  a  differential  element  of  either  the  bourrelet  or 

rotating  band,  it  may  be  written  in  component  form  (relative  to 

S  )  as 
o  ’ 


FL  U 


a 


*  'v,  i- 


(59) 


where  referring  to  Figure  U,  there  results  for  the  bourrelet 

otfv*  =  -  coicc 

ct rkRl(f)l-mf-/~lsi>i«.cosfU<f  J 


and  for  the  rotating  band 

d  £«  “  K  f*i  RJW  COS*.  eLf  1 

d  a  rjRjflL- cos <f  f]~Tsin ¥}<£•? 

-/*jLSiK*cosf]*Tc.os<f}d<f  J 


(61) 


where  and denote  respectively  the  coefficients  of  friction 
at  the  bourrelet/bore  and  ro tating-band/bore  interfaces. 

Hence,  the  resultant  force  applied  to  the  projectile  c.g.  is 


obtained  by  integrating  equation  (59)  over  the  respective 

rotating-band/bore  and  bourrelet/bore  interfacial  contact  surfaces 

defined  by  )  *  <f  *  (f±+tc  )*  1  =  1.  2. 

i  i 

Substituting  equation  (5*0  into  equations  (6o)  and  (6l),  and 
performing  the  required  integrations  for  both  the  rotating  band 
and  bourrelet,  there  results  the  desired  contributions  to  the 
force  components  acting  at  the  projectile  c.g«  (relative  to  Sq) 
due  to  rotating  band  and  bourrelet  interfacial  contact  with  the 
bore,  namely 

f\g--2r;c  o$«.['4l ju.,{  +S,sintft)+  S’  $ -*■<$*  5<  *  )] 

r^^sin^sin^-cos^JfzS^m^S^fc  +  ^n^cos^,)]* 

+  4/[^4sinM.siH^cos^][2S’sin^e4+S2(<fCl*siyi^t^^)b 
-ZT  sin^s/*^} 

-^[^z5mocco<sg+Sii75j][£«S>m^+4^2't-*Si>i^JlC05S?4)]+ 

+2Tccs%3ingA} 

Referring  again  to  Figure  4,  the  moment  arm,  JL-.  ,  from  the 

Mi 

projectile  c.g.  to  the  load  acting  on  the  differential  element  is 
given  by 

iX*  fa *  [ II  Sin^-2.,1-1.  (63) 
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Hence t  the  resultant  moment  applied  to  the  projectile  c.g.  Is 
given  by 


Substituting  equations  (5*0*  (60),  ( 6 1 )  and  (63)  into  equa¬ 
tion  (64),  and  performing  the  required  integrations  for  both  the 
rotating  band  and  bourrelet,  there  results  the  contribution  to 
the  moment  components  acting  at  the  projectile  c.g.  (relative  to 
SQ)  due  to  rotating  band  and  bourrelet  interfacial  contact  with 
the  bore,  namely 


11>  rk  [  W  C.  -  *  rj £J-  ty3-e}  £+  z*, 

Mj. *  coil* »in$h rhjutcos* smf][z //s  in  +Sl(fct+sihficos<£/)l+ 

+2rk/,Tco*%  2^tAneL  ft* 

hi*  4ir*[J,{j*,s*n*3inf-cos1’)+rkf1cos«coiy][2.S'ynft  +  S,(f<  +  5 

♦  )+  co*  5/n^cos^)]4 

*2  rkJttT  sin  H*  sin  V*x  “  *)  taw*  /v,° 


i 


>  (65) 
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The  desired  contribution  to  the  moment  components  M  ,  M  and 

x  y 

M  (written  relative  to  S),  is  obtained  applying  the  transformation 
z 
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mmumm 


'•Hl'T*' 


given  in  equation  (20),  from  which  there  results  noting  Assumption 
(vii) 


m;  -  K/v-M'j.-eM; 


a 


> 


(66) 


m‘,  =  K*  e  m:.  j 

2.3.3  Propellant  Gas  Pressure  and  Compressed 
Air  Loadings 

The  loads  and  moments  applied  to  the  projectile  c.g.  due  to 
propellant  gas  pressure  acting  at  the  base  of  the  projectile  and 
compressed  (’’ram")  air  ahead  of  the  projectile  are  derived  subject 
to  the  following  assumptions: 

(i)  both  the  projectile  base  and  "ram"  air  pressures  are 
assumed  to  be  known  functions  of  time  only,  and  at 
any  instant  are  uniformly  distributed  ovei-  the 
respective  projectile  surfaces  over  which  they  act; 

(ii)  both  the  rotating  band  and  bourrelet  are  assumed  to 
act  as  ideal  seals  (i.e.,  no  gas  leakage); 

(iii)  the  Euler  angles  y>  and  6  satisfy  equation  (3*0. 

Under  Assumptions  (i)  and  (ii),  the  base  and  "ram"  air 

pressures  act  as  effective  hydrostatic  pressures;  generating  a 
force  equal  to  -  p.n«dA  on  the  differential  surface  area  element, 
dA,  with  unit  outward  normal  n.  Further,  from  general  principles 
of  hydrostatics,  these  pressure  loadings  may  be  replaced  by 
resultant  forces  acting  at  the  geometric  center  of,  and  directed 
perpendicular  to,  the  planes  containing  the  bourrelet  and  rotating 
band.  Hence,  the  points  of  application  of  these  resultant  forces 


are  defined  by  the  respective  intoraootion  of  each  of  these 
planes  with  tho  Run  tube  centerline,  as  shown  In  Figure  5* 

Since  *£  Is  parallel  to  the  geometric  axis  of  the  projectile, 
it  follows  that 


(67) 


where  p  and  p^  denote  respectively  the  instantaneoua  "ram"  air 
a  b 

and  projectile  base  pressures,  and  A  denotes  the  projected  area 
of  the  bourrelet  and  rotating  band  planes  perpendicular  to 
Under  Assumption  (iii),  it  follows  that 

Trr^  (68) 


Hence, 


the  resultant  pressure  load, 


is  given  by 


( n-n)rrr^T 

which,  noting  that  in  view  of  Assumption  (iii) 

T  =  t  *  t/.  - 


(69) 


(70) 


may  be  written  in  component  form  (relative  to  Sq)  for  application 
to  equation  (32)  as 


frN  (^-g)Trr/ 
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Projectile  Base  and  "Ram"  Air  Pres 


Once  again  in  view  of  Assumption  (iii),  it  may  be  shown  that 

the  moment  arm,  Jt  ,  from  the  projectile  c.g.  to  the  "ram"  air 

Ql 

load,  F  ,  is  given  by 

a 

(72) 

while  the  corresponding  moment  arm,  ,  to  is 

-  (x^JUe)i 1  (73) 

Hence,  the  resultant  moment  applied  to  the  projectile  c.g. 
is 

=  <7*> 


Substituting  equations  (67),  (72)  and  (73)  into  equation 
(7*0,  and  performing  the  indicated  vector  operations,  there 
results  the  desired  contribution  to  the  moment  components  acting 
at  the  projectile  c.g.  (relative  to  S)  due  to  propellant  gas  and 
"ram"  air  pressures,  in  the  form 


K- 0 

s-irrhll(fi-p)z.,  +  >  (75) 

2.3.4  Summation  of  Applied  Loads  and  Moments 
Summing  corresponding  force  components  from  equations  (43), 
(62)  and  (7l)*  there  results  for  application  to  equation  (32) 
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(76) 


> 

Similarly,  summing  corresponding  moment  components  from 
equations  (66)  and  (75)*  there  results  for  application  to 
equation  (33) 

M,«  fit 

m,=  m:+mi 

2 . 4  Solution  Technique 

Within  the  framework  of  the  assumptions  introduced,  equations 
(32)  and  (33),  with  the  applied  loads  and  moments  as  defined  in 
equations  (76)  and  (77) »  prescribe  projectile  in-bore  motion  in 
terms  of  projectile  design  data,  interior  ballistics  data,  and 
gun  tube  design  and  motion  data.  Of  these  required  data,  gun 
tube  motion  data  are  not  known  a  priori,  and  hence,  equations 
(32)  and  (33)  must  be  solved  simultaneously  with  the  equations 
of  the  gun  dynamics  simulation  code  (References  1  and  2). 

Equations  (32),  (33),  (76)  and  (77)  consist  of  a  system  of 
six  simultaneous  ordinary  differential  equations  for  the  determi¬ 
nation  of  the  six  unknowns  v  *  y  ,  z  ,  'V',  $  .  and  T.  However, 

p  eg  eg  '  ’  ’ 


(77) 
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noting  that  is  to  be  prescribed  by  interior  ballistics  data, 
the  first  equation  in  (32)  may  be  deleted  (or  retained  for  use  as 
a  check)}  rendering  a  system  of  five  simultaneous  equations  for 
incorporation  within  the  gun  dynamics  simulation  code. 

To  incorporate  the  projectile  motion  equations  within  the 
simulation  code,  it  is  noted  that  the  parameters  representing 
gun  tube  motion  in  equations  (32)  and  (33)*  namely,  the  compo¬ 
nents  of  Rq,  their  time  derivatives,  the  Euler  angles,  “Y^*  80  and 
,  and  their  time  derivatives,  are  related  to  gun  tube  motion 
parameters  defined  in  the  simulation  code.  In  addition,  the 
negative  of  the  interfacial  contact  loads  between  the  rotating- 
band  and  bore,  and  bourrelet  and  bore,  as  applied  to  the  projec¬ 
tile,  must  be  applied  to  the  gun  tube.  The  specific  details  of 
incorporating  the  projectile  in-bore  motion  equations  herein 
developed  within  the  gun  dynamics  simulation  code  will  be 
presented  in  detail  in  a  separate  report. 
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3.  COMPARISON  WITH  OTHER  PROJECTILE  DESCRIPTIONS  AND  FORMULATIONS 


The  formulation  herein  developed  is  general  in  that  it 
prescribes  the  full  six  dogree-of-freedom  motion  of  a  projectile 
of  finite  geometry  and  inerti.it  traveling  in  a  flexible  gun  tube. 
The  degrees-of-freedom  selected  correspond  to  three  orthogonal 
translational  motions  of  the  projectile  c.g.  relative  to  the  gun 
tube  axis  and  three  (Eulerian)  rotations  of  the  projectile  about 
its  c.g.  (related  to  projectile  pitch,  yaw  and  roll  motions). 

Since  gun  tube  motion  (which  is  unrestricted  within  the  formula¬ 
tion)  is  not  known  a  priori,  the  projectile  equations  of  motion 
herein  formulated  will  be  solved  simultaneously  with  the  equations 
of  the  gun  dynamics  simulat  Lon  <  ode  previously  developed  (Refer¬ 
ences  1  and  .?). 

A  similar,  but  far  more  restrictive  formulation  has  recently 
been  presented  by  S.H.  Chu  (Reference  3).  Chu  permits  the  projec¬ 
tile  three  degrees-of-freedom;  Iwo  orthogonal  translational 
motions  of  the  projectile  c.g.  relative  to  the  gun  tube  and  one 
rotational  (pitching)  motion.  As  a  consequence  of  neglecting  the 
remaining  degrees-of-frccdom,  the  gun  tube  centerline,  the  projec¬ 
tile  c.g.,  and  the  resultants  of  the  rotating  band  and  bourrelet 
contact-load  distributions  with  the  bore  all  lie  in  the  same 
plane.  Hence,  Chu's  formulation  is  essentially  planar;  whereas, 
the  formulation  herein  presented  is  «>f  general  three-dimensional 
character. 

Another  recent,  but  also  restrictive  projectile  motion 
formulation  has  been  presented  by  H.L,  Langhaar  and  A.P.  Poresi 
(Reference  4).  Langhaar  and  lioresl  present  a  rigorous  kinematical 
description  of  a  point  moving  along  n  time-dependent  spare  run  o , 

Vi 


The  point  is  identified  with  the  gooraietric  center  of*  a  rigid  pro¬ 
jectile.  The  time -dependent  space  curve  is  identified  with  the 
centerline  of  a  flexible  gun  tube.  The  projectile  is  further 
characterized  such  that  its  geometric  center  and  c.g.  coincide 
(which  precludes  the  ability  to  investigate  the  effects  of  mass 
eccentricity),  and  such  that  its  geometric  arts  is  directed  along 
the  instantaneous  tangent  to  the  gun  tube  centerline  (which 
precludes  the  ability  to  investigate  the  effects  of  projectile 
pitch  and  yaw  motions  relative  to  the  gun  tube).  The  projectile 
is  permitted  two  degrees-of-freedom  relative  to  the  gun  tube; 
translational  motion  of  its  c.g.  along  the  gun  tube  centerline 
and  rotational  motion  about  the  centerline  (corresponding  to 
projectile  spin).  Accounting  for  rotary  inertia  of  the  projec¬ 
tile  about  its  spin  axis,  there  results  a  traveling  point-mass 
projectile  load  with  superposed  gyroscopic  couple. 

J.J.  Wu  (Reference  5)  also  adopts  a  traveling  point-mass 
projectile  description,  but  with  superposed  traveling  pitching 
moment  due  to  mass  eccentricty  (while  neglecting  rotary  inertia 
about  the  pitch  axis).  Several  other  investigators  (References 
6  thru  8)  have  adopted  the  simpler  traveling  point-mass  descrip¬ 
tion,  with  and  without  mass  eccentricity  and  projectile  spin 
(while  neglecting  rotary  inertia  about  the  spin  axis). 
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k.  CONCLUSIONS 


A  system  of  simultaneous  differential  equations  has  been 
formulated  which  describes  the  general  motion  of  a  projectile  of 
finite  geometry  and  inei’tia  traveling  in  a  flexible  gun  tube. 

The  formulation  permits  the  projectile  six  degrees-of-freedom 
relative  to  the  gun  tube;  three  orthogonal  translational  motions 
of  the  projectile  c.g.  relative  to  the  gun  tube  axis  and  three 
(Eulerian)  rotational  motions  of  the  projectile  about  its  c.g. 
(related  to  projectile  pitch,  yaw  and  roll  motions).  The  formu¬ 
lation  is  presented  in  terms  of  interior  ballistics  data,  projec¬ 
tile  design  data,  and  gun  tube  design  and  motion  data.  The  formu¬ 
lation  accounts  for  projectile  spin,  mass  eccentricity,  elastic 
deformation  of  the  projectile  rotating  band  and  bourrelet,  and 
pro jectile/bore  interfacial  friction  and  torque  transmission. 
Furthermore,  the  formulation  is  unrestricted  regarding  the  nature 
of  gun  tube  motion. 

The  formulation  herein  contained  has  been  compared  with  pro¬ 
jectile  descriptions  and  formulations  employed  by  other  investi¬ 
gators  and,  based  upon  this  comparison,  is  considered  to  be  the  most 
generally  applicable  formulation  appearing  in  recent  literature. 

Since  gun  tube  motion  is  not  known  a  priori,  and  indeed 
since  projectile  motion  affects  gun  tube  motion  and  vice  versa, 
the  projectile  in-bore  motion  equations  herein  formulated  will  be 
incorporated  within  and  solved  simultaneously  with  the  gun 
dynamics  simulation  code  equations  previously  developed;  permit¬ 
ting  replacement  of  the  traveling  point-mass  projectile  descrip¬ 
tion  presently  contained  within  the  code.  Incorporation  of  the 
projectile  motion  equations  within  the  simulation  code  and 

if; 


simultaneous  solution  with  the  gun  dynamics  equations*  as  well  as 
comparison  of  the  results  obtained  with  previously  obtained 
results  treating  the  projectile  as  a  traveling  point-mass*  will 
be  presented  in  a  separate  report. 
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APPENDIX 

The  instantaneous  angular  orientation  of  the  coordinate 

system  S  ,  with  coordinates  (x  ,y  .z  )  and  unit  triad  (t  ,3\  ,1c  ), 
j  0 »  ‘O  O  O  o  o  o 

•  •  •  I 

relative  to  the  coordinate  system  S  ,  with  coordinates  (x  ,y  ,  z  ) 

.  ^  i  *.  *  . 

and  unit  triad  (i  , j  ,k  ),  is  defined  in  terms  of  the  Euler  angles 
(Vo-Bo’V  depicted  in  Figure  A. 

Referring  to  Figure  A,  the  instantaneous  angular  orientation 
of  Sq  is  achieved  by  subjecting  S  to  the  following  consecutive 
rotations : 

(i)  about  z  ,  bringing  x  to  its  final  elevation, 
and  y  to  its  intermediate  orientation, 

(ii)  about  ,  bringing  £  to  its  final  azimuth,  xq, 

«  / 
and  z  to  its  final  azimuth,  f  ; 

(iii)  <f  ^  about  xq,  bringing  )£  to  its  final  orientation,  y  , 

and  £  to  its  final  elevation,  zq  . 

The  direction  cosines,  J?.  . ,  defining  the  transformation 

^  J 

i 

between  Sq  and  S  are  obtained  by  noting  the  relation  between  the 
unit  triads  depicted  in  Figure  A  subsequent  to  each  consecutive 
rotation,  as  follows: 

(i)  following  the  rotation ^  ,  there  results 

r«  COSTS  $-5in*Z; 

£  =  sinVi  l  *  cosy. 

(it)  following  the  rotation  0q,  there  results 

*  cos $•  t0  •+  $in6* if 
=  -  sind »  l0  +  cos 

A- 1 


(iii)  following  the  rotation  there  results 


y  -  sinii  +  cos£*£ 

Eliminating  the  unit  triad  associated  with  the 

intermediate  coordinates  there  results  the  transforma¬ 

tion  presented  in  equation  (l),  with  the  direction  cosines  as 
defined  in  equation  (2). 

The  relation  between  the  angular  velocity  of  Sq  relative  to 

•  ^ 

S  ,  namely  u )Q>  and  the  Euler  angles,  is  obtained  by  noting  that 

Oo  =  +  d  \  *  fc  t 

and  substituting  the  above  transformations  betweon  unit  triads  to 
obtain  in  the  form  prescribed  in  equation  (  1 1  )  ;  which  renders 
the  components  as  given  in  equation  (l2). 


